Interaction-Induced Spin Polarization in Quantum Dots 



M. C. Rogge 1 ^ E. Rasanen 2 , and R. J. Haug 1 

1 Institut fur Festkorperphysik, Leibniz Universitdt Hannover, Appelstr. 2, 30167 Hannover, Germany 
2 Nanoscience Center, Department of Physics, University of Jyvdskyld, FI-40014 Jyvaskyla, Finland 

(Dated: January 29, 2010) 

The electronic states of lateral many electron quantum dots in high magnetic fields are analyzed in 
terms of energy and spin. In a regime with two Landau levels in the dot, several Coulomb blockade 
peaks are measured. A zig-zag pattern is found as it is known from the Fock-Darwin spectrum. 
However, only data from Landau level show the typical spin-induced bimodality, whereas features 
from Landau level 1 cannot be explained with the Fock-Darwin picture. Instead, by including the 
interaction effects within spin-density-functional theory a good agreement between experiment and 
theory is obtained. The absence of bimodality on Landau level 1 is found to be due to strong spin 
polarization. 

PACS numbers: 73.21. La, 73.23.Hk, 73.63. Kv 



Spin properties of semiconductor quantum dots (QDs) 
are of high interest, as the spin of electrons captured in 
a QD could be used to realize a quantum mechanical 
bit, the core of future quantum information technologies 
p|. In few electron QDs, the electronic spin has been 
successfully implemented, manipulated and read (for a 
review see Q). However, in many electron systems, the 
dynamics are still not well understood, especially in high 
magnetic fields. 

For QDs with many electrons, the simplest theoret- 
ical approximation is done with the so called constant 
interaction (CI) model [!, [3]. This model uses the 
single-electron states, most commonly those of a two- 
dimensional harmonic potential (Fock-Darwin spectrum 
and the many-body effects are included just by a 
constant Coulomb repulsion energy. Despite its simplic- 
ity, many features can be qualitatively explained using 
this model. Among those are the formation of Landau 
levels (that allows to introduce the QD filling factor v), 
the crossing of states leading to zig-zag patterns, alter- 
nating spins and spin flips. Beyond the CI model, self- 
consistent calculations successfully described the regular- 
ity of zig-zag patterns and the electron densities within 
certain Landau levels Q- It was found that Landau 
levels form conductive rings in the dot. Especially for 
4 > v > 2, a central region is formed with Landau level 
1 (LL1) and an outer ring with Landau level (LLO), as 
schematically shown in Fig. [1] 

Further investigations on spin blockade and Kondo ef- 
fect helped to gain knowledge about the spin configura- 
tion in the lowest Landau levels (see, e.g., Refs. |8l-ll3|). 
Calculations by Wensauer et al. [lj] predicted a spin 
polarization in LLO at the v — 2 border as a function of 
the electron number. This was also found experimentally 
[ill [TH . In addition, collective spin polarization on the 
highest occupied Landau level has been found in many- 
electron calculations p]| [H, G3 indirectly supported by 



'Electronic address: roggc@nano.uni-hannover.de 




FIG. 1: (color online) Schematic of a single lateral quantum 
dot connected to two leads (Source and Drain). In a per- 
pendicular magnetic field, Landau levels (LL) appear with a 
specific spin polarization. Landau level is located at the 
edge of the dot, Landau level 1 in the center. Arrows show 
the spin configuration within each Landau level. 

experimental data However, despite all these re- 

sults, the level spectrum of many-electron QDs still leaves 
open questions, and further measurements are required. 

Here we concentrate on the spin configuration of both 
lowest Landau levels in the 4 > v > 2 regime. We 
present conductance measurements on a lateral QD with 
approximately 50 electrons in a perpendicular magnetic 
field. A zig-zag pattern of electronic states is found and 
it is compared with many-body calculations within spin- 
density-functional theory pj] (SDFT). While the states 
in LLO show a striking bimodality that can roughly be ex- 
plained with alternating spin states using the CI model, 
the states of LL1 do not at all follow this model. How- 
ever, the behavior can be understood with the SDFT 
calculations showing different spin configurations in the 
two involved Landau levels as illustrated in Fig. [TJ While 
LLO shows a regular filling of orbitals with alternating 
spin- up and spin-down electrons as expected, LL1 is spin- 
polarized with spin-up electrons only. 

The sample is made on a GaAs/AlGaAs heterostruc- 
ture using local anodic oxidation [I9I - I21I ] and electron 
beam lithography. To allow electronic transport, the QD 
is connected to two leads, source and drain, as schemat- 
ically shown in Fig. [T] A side gate is used to tune the 
QD potential. Details about the sample preparation and 
device properties can be found in Ref. [22j |. 

Figure E] (a) shows a measurement of the magneto- 
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FIG. 2: (color online) (a) Differential conductance G as 
a function of magnetic field B and gate voltage Vg- Six 
Coulomb peaks are visible showing a pronounced zig-zag pat- 
tern, (b) The peak positions are transferred to an energy scale 
with the Coulomb energy removed. A cross pattern is found 
with states going down in energy with increasing field (LLO) 
and states going up in energy (LL1). 



transport of the system. The differential conductance G 
is plotted as a function of the perpendicular magnetic 
field B and gate voltage Vg ■ Six Coulomb peaks are vis- 
ible. The positions of these peaks reflect the chemical 
potentials of ground-state transitions. These positions 
can be roughly understood using the CI model. Accord- 
ing to this model, the potentials include a fixed Coulomb 
repulsion energy due to the electrostatic interaction of 
the electronic charge and a term due to the single par- 
ticle excitation spectrum. The single-particle problem is 
usually calculated describing the confinement of the QD 
as a two-dimensional harmonicpotential, which leads to 
the Fock-Darwin spectrum^, Q. In addition, the elec- 
tronic spin is included via the Zeeman term g*/j,BBs z 
with the gyromagnetic ratio g* , the Bohr- magneton \ib 
and the spin quantum number s z — ±1/2. According 
to this model, the pronounced zig-zag structures appear 
as with increasing magnetic field different states of the 
excitation spectrum are energetically favored. 

The bare excitation spectrum is made visible by re- 
moving the Coulomb energy from the chemical poten- 
tials. This is done in Fig. [5] (b). The peak positions of 
the original Coulomb peaks are transferred to an energy 
scale with the Coulomb energy removed (we did actually 
remove a gate voltage dependent Coulomb energy, as the 
electrostatic electron-electron interaction is reduced with 
increasing Vg due to changes in the size of the dot). Now 
the states from the excitation spectrum can be followed 
over several Coulomb peaks (several electron numbers). 
As a result, a pattern is found with states going up in 
energy (some marked blue) and states going down in en- 
ergy (some marked red). These states can be interpreted 
in terms of Landau levels. States with negative slopes are 
due to transport via LLO, states with positive slopes are 
due to transport via LL1. According to the Fock-Darwin 
model with spin, for each Landau level there should be 
a pairing of every two lines. Two adjacent peaks use 
the same orbital state with opposite spin. Indeed such a 
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FIG. 3: (color online) (a) Measured energetic peak distances 
for states in LLO (triangles) and in LL1 (diamonds). The 
distances for LLO show a bimodal behavior due to spin and 
can be fitted with the Fock-Darwin model (dashed lines, 
fkuo = 1-83 meV, g* — —0.71). The Fock-Darwin spectrum 
using these parameters is shown in (b). Values for LL1 [di- 
amonds in (a)] cannot be explained with the Fock-Darwin 
spectrum. 



pairing is visible for LLO but not for LL1. 

In order to analyze this behavior in detail, for each 
Landau Level the energy distances between adjacent 
peaks are plotted in Fig. [3] (a). The distances for LLO 
(triangles) show the expected bimodal behavior due to 
the electronic spin. The data can be fitted (dashed lines) 
using the Fock-Darwin model with confinement strength 
huiQ = 1.83 meV and g* = —0.71. The corresponding 
Fock-Darwin spectrum using these parameters is shown 
in Fig. [3] (b) around the chemical potential fi4Q (black 
trace). Spin- up and spin-down states are shown as solid 
and dotted lines, respectively However, the large nega- 
tive value for the g*-factor (the sign is given by the peak 
amplitudes due to spin blockade [10J) is a first hint, that 
the description with the CI model is incomplete. Due to 
the confinement potential of the dot, the <?*-factor should 
be more positive than the one of —0.44 for bulk GaAs. 
And finally, the description with the CI model is com- 
pletely insufficient, when data from LL1 are included. 
Within the Fock-Darwin model, the energetic distances 
do not depend on Landau level, but they should be iden- 
tical for LLO and LL1. This is obviously not the case. 
In contrast to the Fock-Darwin spectrum [Fig. |3J (b)] 
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the measured values for LL1 [diamonds in Fig. [3] (a)] do 
not show a bimodal behavior. Moreover, the values are 
larger as for LLO and increase nonlinearly with magnetic 
field. Thus a more sophisticated approach is needed in- 
cluding electron-electron interactions that go beyond the 
CI model. 

In the many-body approach we consider the iV-electron 
Hamiltonian 
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where A is the external vector potential (in symmet- 
ric gauge) of the homogeneous, perpendicular magnetic 
field B = Bz, and the last two terms correspond to 
the external potential in the harmonic approximation 
14 xt (r) = rrCuJ^r 2 jl with hjj a — 4 meV, and the Zee- 
man energy E z = g* hbBs z , respectively. We apply here 
the conventional effective-mass approximation with the 
GaAs material parameters: m* = 0.067 and e* = 12.4. 
For the gyromagnetic ratio we have chosen g* = —0.30. 

We solve the ground-state energies associated with the 
iV-electron Hamiltonian (pj by applying the SDFT [l8[ 
in the collinear-spin representation. We note that the 
external vector potential is retained in the corresponding 
Kohn-Sham Hamiltonian. To approximate the exchange- 
correlation energy E xc we use the local spin-density ap- 
proximation (LSDA) with a paramctrization of the corre- 
lation energy in the homogeneous 2D electron gas by At- 
taccalite et al. [23| For total-energy calculations on QDs, 
SDFT with the LSDA has been shown to be a reliable 
scheme in comparison with quantum Monte Carlo cal- 
culations, even in relatively high magnetic fields [l3|, [13, 
[24|. Furthermore, it has been found that the current- 
SDFT [25] (with the local-vorticity approximation) does 
not lead to a considerable improvement over the SDFT 
results [24]. In the numerical calculations we apply the 
octopus code package [26j [. 

Figure LU[a) shows the chemical potentials /ijv = -EW — 
E N _i of the SDFT calculations for TV = 46 ... 48. The 
resulting energy-level distances extracted from the chem- 
ical potentials are shown in Fig. |U[b). For LL1 we find a 
clear single mode, whereas LLO has a bimodal behavior as 
a function of B. This is in a good qualitative agreement 
with the experimental result in Fig. Ufa). 

To analyze the physics behind the energy-level spac- 
ings, we focus in the following on the spin configurations 
given by the SDFT calculations. Their connection to the 
chemical potentials is shown in Fig. [5] A typical por- 
tion of a Coulomb peak in the zig-zag regime is shown 
together with the spin configurations for both Landau 
levels for the adjacent regions of constant electron num- 
ber. LLO shows alternating spin filling resulting in the 
bimodality of peak sections with negative slopes. The 
bimodality is due to the energy difference in adding an 
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FIG. 4: (color online) Results of the SDFT calculations: (a) 
chemical potentials /M6, ^47 and /i48 without Coulomb energy. 
As in the experiment, a cross pattern is found, (b) Energetic 
distances for the two Landau levels. The data qualitatively fit 
to the measurements. For LLO a bimodal behavior is found 
(squares). The values for LL1 (diamonds) increase with in- 
creasing field and do not show a bimodality. 
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FIG. 5: (color online) Spin configuration resulting from the 
SDFT calculations for a typical section of the zig-zag pattern. 
While electrons in LLO are filled in regularly with alternating 
spin, LL1 is spin polarized with electrons having spin up only. 



electron with spin up to a new orbital compared with 
adding an electron with spin down to a half-filled or- 
bital. We note that the energy difference decreases as a 
function of B. This is due to both the Zeeman effect and 
the fact that the increasing magnetic confinement brings 
the energy states in LLO closer to each other, so that the 
energy in occupying a new orbital decreases with respect 
to the filling of a half-filled orbital (costing Coulomb en- 
ergy). 
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In contrast with the situation for LLO, the higher Lan- 
dau level (LL1), corresponding to peak sections with pos- 
itive slopes in Figs. 0] and EJ is completely spin-polarized 
with spin up. Thus, the filling mechanism is always the 
same and no bimodality is found. The single mode goes 
up in energy due to the effective increase in the level 
spacings as a function of B, i.e., LL1 becomes more com- 
pact, so that the filling of the spin- up orbitals becomes 
energetically more costly. 

We point out that the qualitative features in Figs. 0] 
and [5] arc stable with respect to the values chosen for the 
gyromagnetic ratio. The reason behind the strong spin 
polarization of LL1 is the high density of states close to 
the Fermi energy leading to collective, local ferromag- 
nctism familiar from the Stoncr effect and Hund's 
rule. Energetically, this many-body transition corre- 
sponds to gaining electronic exchange at the expense of 
kinetic energy in occupying higher levels. Similar be- 
havior corresponding to the formation of "spin-droplets" 
has been predicted by quantum Monte Carlo and SDFT 
calculations on QDs, indirectly supported by the experi- 
mental data of Coulomb-blockade peak positions [l3[ . It 



was concluded that the spin-droplet formation requires 
a relatively large number of electrons (N > 30) and 
considerable strength of electron-electron interactions. 
These conditions are fulfilled by our device, and hence 
the Landau-level distances visible in our excitation spec- 
tra can be considered as a direct and distinct evidence 
of the spin polarization on the highest occupied Landau 
level. 

To conclude, we analyzed the spin configuration of a 
single quantum dot containing approximately 50 elec- 
trons. In high magnetic fields electronic states were in- 
vestigated for the two lowest Landau levels. For LLO dis- 
tances in the excitation spectra show a bimodal behavior, 
while for LL1 only one mode appears. This is explained 
with two different spin configurations in the two Landau 
levels. While there is no spin polarization in LLO, LL1 is 
completely spin polarized. This is a direct confirmation 
of the interaction-induced collective spin polarization in 
quantum dots. 

This work has been supported by BMBF via 
nanoQUIT, by QUEST, and by the Academy of Finland. 
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